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Abstract: Super symmetric extensions of the Standard Model that incorporate the 
axion solution to the strong CP problem necessarily contain also the axino, the 
fermionic partner of the axion. In contrast to the neutralino and the gravitino, the 
axino mass is generically not of the order of the supersymmetry-breaking scale and 
can be much smaller. The axino is therefore an intriguing candidate for a stable su- 
perpartner. In a previous Letter [1] it was shown that axinos are a natural candidate 
for cold dark matter in the Universe when they are generated non-thermally through 
out-of-equilibrium neutralino decays. Here, we extend the study of non-thermal pro- 
duction and include a competing thermal production mechanism through scatterings 
and decays of particles in the plasma. We identify axino masses in the range of tens 
of MeV to several GeV (depending on the scenario) as corresponding to cold axino 
relics if the reheating temperature Tr is less than about 5 x 10 4 GeV. At higher Tr and 
lower mass, axinos could constitute warm dark matter. In the scenario with axinos 
as relics, the gravitino problem finds a natural solution. The lightest superpartner of 
the Standard Model spectrum remains effectively stable in high-energy detectors but 
may be either neutral or charged. The usual highly restrictive constraint Q x h 2 < 1 
on the relic abundance of the lightest neutralino becomes void. 
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1. Introduction 

Among several possible ways of addressing theoretical puzzles of the Standard Model, 
two are generally believed, for quite independent reasons, to be very compelling. 
Softly-broken low-energy supersymmetry (SUSY) seems to be the most attractive 
way of solving the hierarchy problem and of linking the electroweak scale with physics 
around the Planck scale. Moreover, it has the generic property of contributing only 
small corrections to electroweak parameters and also of predicting a light Higgs 
boson, both in agreement with the outcome of current precision measurements. 

The most compelling way of resolving the strong CP problem of QCD seems 
to be provided by invoking the Peccei and Quinn (PQ) mechanism [2]. There, the 
CP- violating 6*-term in the QCD Lagrangian, which is experimentally required to be 
excessively tiny, is replaced by a term involving a new fundamental pseudo-scalar 
field, the axion. This is achieved by introducing a global, chiral U(l) symmetry 
group, which is spontaneously broken at some high energy scale /„ ~ 10 11 GeV. The 
QCD anomaly breaks this chiral U(l) symmetry at the one-loop level, and hence 
the axion becomes not a true Goldstone boson but a pseudo-Goldstone boson with 
a tiny mass of order AQ CD // a [3]. As the lightest pseudo-scalar consistent with the 
known particle phenomenology, the axion can be very important in cosmology and 
astrophysical processes. Axion physics and cosmology are very rich fields which have 
been thoroughly explored in the literature [4] . 
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Axions can also be the dark matter (DM) in the Universe. Indeed, being ex- 
tremely light and feebly interacting, and having a lifetime of order > 10 40 years, 
axions can be considered as stable for all practical purposes. Therefore, if axions are 
copiously produced in the early Universe, they can contribute substantially to the 
relic mass density. Since they are spin-zero particles, we can consider hot axions and 
cold axions separately. Hot axions are produced when the PQ symmetry is broken at 
T < f a - Such hot axions are diluted by inflation if the reheating temperature after 
inflation is less than about 10 9 GeV. At low temperature it is difficult to produce 
hot axions because of their tiny interaction strength. Since the axion potential is ex- 
tremely flat, the axion vacuum expectation value (VEV) (a) will not move until the 
temperature of the Universe falls below ~ 1 GeV. These coherent axions produced 
around the QCD phase transition have very small kinetic energy and are therefore 
cold. They can constitute the DM of the Universe if /„ ~ 10 12 GeV [5]. 

An inevitable prediction of combining the two well-motivated and independent 
hypotheses of axions and supersymmetry is the existence of the SUSY partner of 
the axion, the axino. Being massive and both electrically and color neutral, axinos 
are an intriguing possibility for the lightest supersymmetric particle (LSP) and a 
WIMP. In the presence of i?-parity conservation there may be dramatic cosmological 
implications of the existence of such stable massive particles. In particular, they 
may contribute substantially to the relic mass density and may constitute the main 
component of the DM in the Universe. This is the possibility that we will explore in 
this study. 

There are two other SUSY particles with well-known, and much-studied, cosmo- 
logical properties which, ever since the early days of low-energy SUSY, have been 
considered for DM candidates. Perhaps the most popular of them is the lightest 
neutralino, Xi which, if stable, often provides the desired amount of relic abundance 
for natural ranges of other superpartner masses. Furthermore, since its characteristic 
interaction strengths are often of a sizeable fraction of the electroweak interactions, 
the neutralino has a good chance of being detected in high energy colliders and, if it 
indeed constitutes the dark matter, in WIMP dark matter searches. 

The other traditional candidate, the gravitino, G, arises by coupling SUSY to 
gravity in supergravity or superstring scenarios. It has been long known that grav- 
itinos could be copiously produced in thermal processes in the early Universe and 
often suffer from the well-known "gravitino problem" [7] of excessive destruction of 
light elements and of contributing too much to the energy density. These prob- 
lems can be avoided if the reheating temperature after inflation is sufficiently low, 
Tr < 10 9 GeV [6, 7, 8]. Therefore, the gravitino remains a possibility for dark matter. 
(For recent work see Refs. [9, 10, 11].) 

In contrast to the neutralino and gravitino cases, axinos have, undeservedly, 
attracted much less attention in the literature. This is even more surprising in view 
of the fact that axinos possess some distinctive properties that lead to important 
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cosmological implications. In particular, in contrast to the case of the gravitino as 
well as the neutralino and other ordinary superpartners, the mass of the axino is 
generically not of the order of the SUSY-breaking scale and, as we will see later, can 
be much lower [12, 13, 14]. This alone makes the axino an intriguing candidate for 
the LSP and dark matter. 

One important exception was provided by a comprehensive study of Rajagopal, 
Turner and Wilczek [13]. These authors considered both light (~ keV range) and 
more massive (~ GeV) axinos and studied their cosmological implications as either 
the LSP or the next-to-LSP (NLSP). They derived an upper bound < 2keV on 
the mass of stable, "primordial" axinos which would hold in the absence of inflation. 
They also concluded that, in a class of interesting models, the axino could have a 
mass satisfying this bound and would constitute so-called 'warm ' dark matter. Such 
light axinos were studied also by other early papers [15, 16]. 

While warm DM has some interesting properties (see, e.g., Ref. [17]), it has been 
the standard cosmological lore of the last several years that the invisible dynamical 
component of the mass-energy density of the Universe is probably predominantly in 
the form of cold dark matter (CDM) [18]. The paradigm has been based on good 
agreement of numerical simulations of the formation of large scale structures with 
observations [18]. More recent studies have revealed some possible problems with 
predictions of the standard CDM theory on sub-galactic scales, which still need to 
be clarified. Nevertheless, we believe that it is still well-justified to assume that most 
of the DM in the Universe is cold. 

In a previous paper [1], we pointed out that axinos can naturally form cold dark 
matter. One way of producing such CDM axinos is through out-of-equilibrium de- 
cays of a heavy enough superpartner. Because axino couplings to other particles 
are suppressed by l/f a , as the Universe cools down, all heavier SUSY partners first 
cascade-decay to the lightest ordinary SUSY particle (LOSP), e.g., the lightest neu- 
tralino. (By "ordinary" we mean any of the superpartners of the Standard Model 
particles.) The LOSPs then freeze out of thermal equilibrium and subsequently de- 
cay into axinos. If LOSP is a neutralino of tens of GeV or heavier, its lifetime, which 
scales like l/m x , is often much shorter than 1 sec and the decay takes place before 
Big Bang Nucleosynthesis (BBN) [1]. 

It is worth mentioning that the neutralino will nevertheless appear to be stable 
in high-energy colliders since its lifetime will typically be significantly larger than 
~ 10~ 7 sec. The same would be true if the LOSP carried electric charge. In this case 
the LOSP would appear in a detector as a stable, massive, charged state [13]. 

Furthermore, since after freeze-out all neutralinos convert into axinos, one ob- 
tains a simple expression for the axino relic abundance Vtah 2 = (ma/m x ) Q x h 2 [1]. 
(Here Qa = ps/Pcrit, where is the relic density of axinos and p cr i t is the critical 
density, h is the dimensionless Hubble parameter, mj and m x are the respective 
masses of the axino and the neutralino, and Q x h 2 is the abundance the neutralinos 
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would have had today.) It follows that cases where Quh 2 ~ 1 will typically corre- 
spond to supersymmetric configurations for which Q x h 2 ^> 1. In other words, cases 
traditionally thought to be excluded by the condition Q x h 2 < 1 will now be readily 
allowed (actually, even favored!). They would also typically correspond to larger su- 
perpartner masses, which may be unwelcome news fo collider searches, but appears 
to be favored by gauge coupling unification [19, 20] as well as flavor and proton decay 
constraints [21, 20]. 

In addition to the non-thermal production (NTP) discussed above, axinos can 
also be generated through thermal production (TP), namely via two-body scatter- 
ing and decay processes of ordinary particles still in thermal-bath. (We stress that, 
despite the name of the process, the resulting axinos will typically be already out of 
thermal equilibrium because of their exceedingly tiny couplings to ordinary matter, 
except at very large Tr. They will nevertheless be produced in kinetic equilibrium 
with the thermal-bath and thus their momenta will have a thermal spectrum inher- 
ited from the scattering particles in the plasma.) While the analogous process of 
gravitino production has been studied in the literature [8, 9], this has not yet been 
done in the case of axinos. 

In this study we aim to provide a detailed exploration of both TP and NTP 
mechanisms of producing axinos through processes involving ordinary superpartners, 
and to compare their relative effectiveness. Both mechanisms are quite generic, but 
specific results will be somewhat model-dependent. For definiteness, we will be 
working in the context of the Minimal Supersymmetric Standard Model (MSSM). 
Also, while any ordinary superpartner (neutral or not) could in principle be the 
LOSP, we will concentrate on the case of the neutralino. We will also, for reference, 
often assume f a = 10 11 GeV, although it will be relatively straightforward to rescale 
our results for other values of f a . Furthermore, since the axions and axinos are 
produced through different mechanisms and both can be CDM, one could consider 
the possibility that both contribute sizeably to the CDM relic abundance. While 
this scenario remains a viable possibility, we will, for definitness, explicitly assume 
that it is the axino that dominates the Universe. 

We will show that, unsurprisingly, TP is more important at larger values of the 
reheating temperature T R . On the other hand, NTP will dominate at rather low (but 
not unreasonably low) values T R < 5 x 10 4 GeV. This is also, broadly, the range of 
Tr for which other constraints will be satisfied. 

Clearly we are most interested in the cases where the relic density of axinos will 
be close to the critical density p CT it- More precisely, we will require 0.1 < Qah 2 < 0.3. 
This condition will provide a strong constraint on the scenario. In particular, in 
the NTP case of axino production through ~ (9(100 GeV) bino decay, it will imply 
rria > C(10MeV) to 0(1 GeV), depending on whether we allow the SUSY-breaking 
scale M SUSY to take very large values (tens of TeV) or require M SU sy ^, 1 TeV, 
respectively. 
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Axinos will initially be nearly always relativistic (unless they are nearly degen- 
erate with the parent neutralinos) , hence they have to be heavy enough to "cool 
down" and become non-relativistic, or cold, by the epoch of matter dominance. We 
find that this leads to a somewhat weaker condition m a > (9(100 keV). 

A similar bound arises from requiring that the relativistic axinos do not con- 
tribute too much to the energy density of radiation during BBN, since their pro- 
duction will take place near the time of BBN. Furthermore, a significant fraction 
of neutralino decays will proceed into gg-pairs (+axino) through an intermediate 
photon or Z-exchange. The resulting hadronic showers may cause excessive destruc- 
tion of light elements. We find that this can be easily avoided but that the specific 
bounds are strongly model-dependent. It would suffice to simply assume a large 
enough neutralino mass m x for the decay to take place before BBN. We find that 
such a crude requirement would be unnecessarily over-constraining. Nevertheless, for 
smaller m x the resulting lower bound on m„ will indeed be typically much stronger 
than the above bounds. For example, for a bino of 60 GeV one finds m„ > 360 MeV. 
However, increasing m x to 150 GeV removes the bound altogether. 

We should mention other non-thermal mechanisms for re-populating the Universe 
with axinos and/or other relics such as gravitinos and moduli. First, it has recently 
been claimed [22] that gravitationally interacting particles can be copiously produced 
through inflaton field oscillations in the reheating or preheating process, and that 
non-thermal production of gravitinos can be much more efficient than through ther- 
mal processes. 1 This would further aggravate the gravitino problems with nucleosyn- 
thesis, thus leading to a much more severe constraint on the reheating temperature, 
sometimes as low as T R < 10 5 GeV, depending on the model. We do not expect 
the mechanism to contribute to the production of axinos in any significant way, es- 
pecially in the regime of low T R mentioned above. Furthermore, such mechanism 
would strongly depend on the inflationary model and on the interactions between 
the inflaton, and the axino and it would seem difficult, if not impossible, to obtain 
the required abundance of DM. 

Second, axinos can be produced in gravitino decays, along with their non-SUSY 
partner, the axion [25]. It is worth emphasizing that the gravitino problem can 
be easily resolved if the axino is the LSP and the hierarchy of masses ttilosp > 
mg > rria is assumed. The process will take place long after nucleosynthesis (t ~ 
10 8 sec) but decay products will now be completely harmless. The lightest ordinary 
superpartner will also decay directly to axinos, thus by-passing the dangerous late 
decays to gravitinos and energetic photons. In the relatively low Tr < 10 5 GeV 
regime, where NTP of axinos dominates, there is not even the gravitino problem 
anymore [22]. However, the mechanism resolves the gravitino problem even at much 
larger T R < 10 15 GeV studied in Ref. [25]. In other words, it seems that the gravitino 

1 This claim has very recently been disputed [23]. See also Ref. [24]. 
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problem can be resolved, assuming the above hierarchy, in a model-independent way. 
At large Tr the axino would not contribute to cold DM. Instead, it would be warm 
DM while cold DM could be provided by, for instance, the axion. 

It is clear that in discussing cosmological implications of axinos one quantity of 
crucial importance is the axino mass itself. Unfortunately, unlike the axino coupling, 
this quantity is rather poorly determined and is strongly model-dependent. As men- 
tioned above, a distinctive feature of axinos is that typically their mass is not set 
by the SUSY-breaking scale and therefore can be much lower. In fact, it can span 
a wide range from ~ eV to ~ GeV [14], depending on the model. In this study we 
will therefore treat basically free parameter. 

Finally, we briefly mention the saxion - the R = 1 spin-zero component of the 
axion supermultiplet. 2 As with other scalars, the mass of the saxion arises from a 
soft term and is typically of the order of the SUSY-breaking scale. Cosmological 
properties of saxions can be quite important [26, 27]. In particular, saxions decay 
relatively fast which may cause significant entropy production which would reheat 
the Universe and lead to other consequences. However, these effects will not be 
relevant to us if the saxion mass is below 1 TeV which we will assume here. We will 
occasionally comment on the saxion below when relevant for axino cosmology. 

The plan of the paper is as follows. In Section 2 we discuss relevant axion and 
axino properties. In particular, we concentrate on the axino mass and couplings. In 
Section 3 the thermal production of axinos is analyzed in detail. Section 4 summarizes 
and extends the previously obtained results in the case of NTP [1]. Astrophysical and 
cosmological constraints are discussed in Section 5. In Section 6 the two production 
mechanisms are compared and further discussed. Finally, in Section 7, implications 
for cosmology and for collider phenomenology are briefly discussed and concluding 
remarks are made. 



2. The axino 

First, let us briefly summarize the main properties of axions. The axion can be 
defined as a pseudo-scalar particle with the effective interaction given by 

where M is a model-dependent mass scale, F is the field strength of the gluon field 
and F is its dual. The potential arising from the above interaction settles (a) to zero. 
In general, there can arise some other small terms in the axion potential in addition 
to the one arising from (2.1). (When CP violation is considered (a) = is shifted, 

2 In the literature s is called saxion (= scalar partner of axion) [13], or saxino (= the scalar partner 
of axino) [27]. Since superpartner names ending with "-ino" seem to be reserved for (Majorana) 
fermions, here we will use the name 'saxion'. 
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but the shift is extremely small since a square of the weak interaction coupling and 
the pion mass appear in the resulting expression [28].) At low energy, one may not 
question how M is generated. It can arise in renormalizable or non-renormalizable 
theories, or in composite models. If a fundamental theory such as superstrings allows 
the above non-renormalizable interaction, then M is the compactification scale [29]. 
If a confining force at a scale A a produces the Goldstone boson with the above 
coupling at low energy, then M is of order A a and the axion is composite [30] . If the 
axion resides in the complex Higgs multiplet(s), it can be derived through the PQ 
mechanism with the spontaneously broken chiral U(l) global symmetry. 

In phenomenologically acceptable axion models, the mass scale M is given by 
M = Mf a /a s , where f a is the PQ scale, 10 9 GeV < f a < 10 12 GeV. In these models 
the axion mass is given by m a ~ ^qcd/ fa ~ lCT 2 -lCT 5 eV; in other words the axion 
is very light. These models can be classified as hadronic axion models (usually called 
the KSVZ models [31]) and DFSZ axion models [32]. Let us briefly recall some 
basic features of these two most popular implementations of the PQ mechanism. 
(Detailed presentations can be found in several excellent reviews, e.g. in [4].) In 
both one assumes the existence of a complex scalar field 0, which is a singlet under 
the SM gauge group but carries a PQ charge. When develops a VEV, ((f)) = f a , 
the PQ symmetry gets broken and the phase of becomes the axion field. To relate 
the global charge of the complex scalar to the PQ charge (i.e. to couple the axion 
to the gluon anomaly), in the KSVZ scheme [31] one further introduces at least one 
heavy quark Q 3 which couples to through the term 

Uq = Iq<PQrQl + h.c. (2.2) 

where Jq is the Yukawa coupling. The PQ charge of is +2, while those of Qr^l 
are ±1. The breaking of the PQ symmetry gives large mass to the heavy quark, 
m Q = /q/o- The axion interacts with ordinary matter through loops involving the 
exchange of the heavy quark (i.e. through the anomaly term). 

This leads, below the PQ scale, to effective axion interactions terms; among these 
the most important one is its interaction with gluons 

£a 9g = ^FF. (2.3) 

In the DFSZ model [32], instead of the heavy quark, one introduces two Higgs 
doublets [2, 33] Mj d of SU(2)l x U(1)y which couple the SM quark sector to 0. 
The doublets carry hypercharges ±1/2 and PQ charges —Q u ,d, respectively. The n 
field couples only to up-type quarks, while 0^ couples only to down-type quarks and 
leptons which carry PQ charges so as to respect the PQ symmetry. 

3 In general one can introduce, instead of a single colored fermion, a multiplet, e.g. a non-singlet 
vector-like quark representation of SU(2) x U(l)y 
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After the breaking of the PQ symmetry, again effective axion interactions with 
ordinary matter are generated, including the term (2.3). In both models there ap- 
pears an additional quantity, the number of different vacua (or, equivalently, domain 
walls) N: N = 1(6) for the KSVZ (DFSZ) model. For our purpose it is enough to 
note that its effect will be to replace f a by f a /N. 

Both models can be readily supersymmetrized [34, 35]. The scalar field be- 
comes promoted to the superfield, and accordingly the other fields as well. In par- 
ticular, the axion, being the phase of 0, is also promoted to an axion supermultiplet 
$ consisting of the pseudo-scalar axion a, its fermionic partner, the axino a, and the 
scalar partner, the saxion s 

$ = (s + ia) + \fTaQ + F*66. (2.4) 
v2 

The axino is thus a neutral, R — —1, Majorana chiral fermion. Adding supersymme- 
try opens up a wide range of choices for the implementation and breaking of the PQ 
symmetry through a choice of different superpotentials and SUSY-breaking schemes, 
as will be illustrated below. 

Now let us move on to the discussion of axino masses. We will briefly summarize 
the relevant results of several previous studies [12, 15, 13, 36, 14, 37]. The overall 
conclusion will be that the mass of the axino is strongly model-dependent [13, 14]. 
It can be very small (~ eV), or large (~ GeV), depending on the model. 

It is worth emphasizing that, unlike the case of the gravitino and ordinary su- 
perpartners, the axino mass does not have to be of the order of the SUSY-breaking 
scale in the visible sector, M SUSY [12, 13, 14]. In global SUSY models, one sets 
-Msusy = 0(1 TeV) on the basis of naturalness. In SUGRA models M SU sy is set 
by the gravitino mass mg ~ Mg/Mp, where Ms ~ 10 11 GeV is the scale of local 
SUSY-breaking in the hidden sector and Mp ~ 2.4 x 10 18 GeV denotes the Planck 
mass. 

It is easy to see why mj is not generically of the order of M SU sy- In the case 
of unbroken SUSY, all members of the axion supermultiplet remain degenerate and 
equal to the tiny mass of the axion given by the QCD anomaly. Once SUSY is 
broken, the saxion, being a scalar, receives a soft-mass term m 2 s s 2 [12, 38], where 
m s ~ Msusy, similarly to the other scalar superpartners. 

Of course, the axion, being instead a phase of the fields whose VEVs break 
the PQ symmetry, does not receive soft mass. Likewise, one cannot write a soft 
mass term for the axino since it is a chiral fermion (for the same reason there are 
no soft terms for, for instance, the MSSM higgsinos) and a superpotential term 
W ~ (mass parameter) • <3><3> is absent due to the PQ symmetry. The lowest-order 
term one can write will be a non-renormalizable term of dimension-5. The axino 
mass will then be of order M| USY // a ~ IkeV [12]. (A more detailed, and more 
sophisticated, explanation can be found, e.g., in Ref. [14].) 
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However, in addition to this inevitable source of axino mass, one can easily 
generate much larger contributions to at one-loop or even at tree-level. One-loop 
terms will always contribute but will typically be < Msusy (KSVZ) or even <C Msusy 
(DFSZ). Furthermore, in non- minimal models where the axino mass eigenstate comes 
from more than one superfield, arises even at the tree-level. In this case can 
be of order M SUSY but can also be much smaller. These points are illustrated by the 
following examples. 

First let us see how tree-level contributions can be generated from superpotentials 
involving singlet fields. For example, consider the case where the PQ symmetry is 
assumed to be broken by the renormalizable superpotential term (in KSVZ or DFSZ 
models) [39] 

W = fZ(S 1 S 2 -f 2 ), (2.5) 

where / is a coupling, and Z, Si and S 2 are chiral superfields with PQ charges of 0, +1 
and —1, respectively. In this case the axino mass can be at the soft SUSY-breaking 
mass scale [36]. It arises from diagonalizing the mass matrix of the fermionic partners 
Si, S 2 , and Z, 

j m~ a ff a \ 

\m~a ff a \, (2.6) 

\f Jaffa / 

where = f(Z) and ff a ~ 10 11 GeV. The mass matrix in Eq. (2.6) is for three two- 
component neutral fermions. These three will split into one Dirac fermion, and one 
Majorana fermion which is interpreted as the axino. The corresponding eigenvalues 
are A = -m~ a and A = ±y/2ff a + 0(m~ a ). In the global SUSY limit, (Z) = and the 
tree-level axino mass would be zero. However, when Si and S2 acquire VEVs and 
soft terms are included, V = |/| 2 (|Si| 2 + \S 2 \ 2 )\Z\ 2 + (AJS^Z - A 2 f f 2 a Z + h.c), 
a linear term in Z is generated which induces (Z) of order (A 1 — A 2 )/ f where A i2 
are the soft trilinear mass parameters. The axino mass thus arises at the soft mass 
scale [36]. 

However, by choosing a more complicated superpotential, one can significantly 
lower the axino mass [14] even at the tree-level. Consider a supergravity superpo- 
tential consistent with the PQ symmetry as 

W = fZ(S 1 S 2 - X 2 ) + ^A(X - M) 3 , (2.7) 

where X carries a zero PQ charge. In this case a minimization of the potential 
resulting from W is much more complicated, and an approximate solution gives 
the lightest eigenvalue of the fermion mass matrix with (V) = which is m a = 
0(A — 2B + C) + 0(rriy 2 / f a ), where A, B and C are the respective trilinear, bilinear 
and linear soft breaking parameters. For the standard pattern of soft breaking terms, 
B = A — m 3 / 2 and C = A — 2m 3 / 2 , the leading contribution A — 2B + C vanishes 
and the tree-level axino mass becomes of order mL/ '/„ ~ 1 keV. 
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The above example shows that in calculating even the tree-level axino mass one 
must carefully consider the full potential V generated in supergravity models. In 
general, the tree-level axino mass either can be of order 7713/2 or, depending on the 
PQ sector of the model and the pattern of soft breaking parameters, can be much 
less, as shown in the above examples. The detailed conditions for this to happen 
were analyzed by Chun and Lukas [37]. In particular one can also have a more 
complicated PQ sector and break the PQ symmetry by non-renormalizable terms in 
the superpotential leading to interesting cosmological implications [40]. 

If the tree-level axino mass is either zero or of order mL/ f a , the contribu- 
tion from loop diagrams can become more important. In the global SUSY ver- 
sion of the KSVZ model, axino mass will arise at the one-loop level with a SUSY- 
breaking A-term insertion at the intermediate heavy squark line. Then one finds 
m a ~ (/q/8tt 2 )^4, where fg is the Yukawa coupling of the heavy quark to a singlet 
field containing the axion (compare Eq. (2.2)), which gives m 5 in the range of < a 
few tens of GeV [41, 36]. In the DFSZ model [32] (in either global SUSY or super- 
gravity), where such A-term contribution is absent, the axino mass remains in the 
keV range, as was pointed out in [15, 13]. Thus, in some models the axino mass can 
be rather small. 

In gauge-mediated SUSY-breaking models (GMSB), the pattern of the axino 
mass is completely different. In the GMSB approach, the tree-level effect of super- 
gravity is not important since the gravitino mass is much smaller than soft breaking 
masses, typically of order ~ eV to ~ keV. These models attracted much interest 
in the recent past since the gravitino production rate is quite large which leads to 
important implications for cosmology [42] and accelerator experiments [43]. In the 
GMSB models, the axino mass is further reduced by one more power of f a and can 
range from 10~ 9 eV to 1 keV [44]. It remains model-dependent but again can be 
smaller than the gravitino mass. We will not discuss GMSB cases here any more. 

In conclusion, we see that a complete knowledge of the superpotential and 
supersymmetry-breaking mechanism is necessary to pin down the axino mass. In 
general it can range from ~ eV to ~ tens of GeV. For the sake of generality, in 
discussing cosmological properties of axinos, we will use as a free parameter. 

In addition to its mass, the axino couplings to gauge and matter fields are of 
crucial importance for the study of its cosmological abundance. In general, the 
couplings of the members of the axion/axino supermultiplet will be determined by 
the PQ symmetry and supersymmetry. The most important term for our purpose 
will be the coupling to gauge multiplets. It is given by the same interaction term as 
the one which gives rise to the "0" term in the QCD Lagrangian. In supersymmetric 
notation it can be written as 



L = 



OtyCaYY 



[B -> W]l 




[B->F]+h.c, 
(2.8) 



4V2n(f a /N) 



AV27i(f a /N) 
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where ai {% — Y, 2, s) are the coupling constant strengths of the SM subgroups, and 
B, W3 and F are the respective gauge supermultiplets. The coefficient C agg = 1 
is universal, while C a yy and C a ww are model-dependent. For processes that in- 
volve the electroweak gauge bosons and superpartners, it is always possible to rotate 
away the interaction of the axion with the SU(2)l gauge multiplet through some 
anomalous global chiral U(l) x rotation since the PQ symmetry is global. Then the 
only interaction that remains to be specified will be the one between the axino and 
the U(l)y gauge multiplet. This will be equivalent to assigning PQ charges to left- 
handed doublets in such a way as to obtain C a ww — 0. One will then also find 
C a YY — C ail . 

The rationale for removing the coupling C a ww is as follows. We would like 
to keep the gauge couplings and top quark Yukawa coupling but neglect the other 
Yukawa couplings. This is because the neglected Yukawa couplings will not be impor- 
tant in the production cross section. Moreover we expect that the lightest neutralino 
is likely to contain a significant bino component [45, 46, 19], so that phenomenolog- 
ically the interaction between the U(1)y field strength and the axion multiplet will 
be important. Keeping this in mind, we will work with the simplest interactions and 
use the basis where the C a ww coupling is absent. 

The procedure of global U(l)x rotation and charge redefinition is somewhat 
subtle and requires a word of clarification on which we digress here. This rotation 
applies only to the high energies of interest to us, but not to the low-energy phe- 
nomena below the QCD scale. We can use a global U(l) x rotation to redefine the 
PQ charges of the quarks like Q' = Q — Q(quark doublet). This will cancel the PQ 
charge of the quark doublets, but will also shift the PQ charges of the other fields. 
Since the Q' charges of the quark doublets are absent, one would expect that there 
is no C a ww coupling either. However, there will be (in DFSZ models) lepton and 
Higgs doublets whose Q' charges will in general be non-zero. Thus the Q' charge 
should instead be defined by the rotation removing C a wwi i n which case Q' charges 
of quark doublets and lepton and Higgs doublets may in general be different from 
zero. Because the top quark Yukawa coupling is kept, the global U(l)x transforma- 
tion must either leave the quark fields unchanged or operate on both the left-handed 
and right-handed quarks and the Higgs field, similarly to the Standard Model U(1) Y - 

The rotation will not change the C agg coupling since U(l)y is anomaly-free. 
On the other hand, for leptons the Yukawa coupling can be neglected and we can 
therefore use any global U(l)x such that the final C a ww coupling is absent. The 
above rotation amounts to changing the PQ charges of the lepton doublets and hence 
to modifying all the couplings of the a- field with the field strengths of and B^. 
If we could neglect also the top quark Yukawa coupling, we could consider instead a 
basis without the C agg coupling where the interaction (2.1) would be absent and no 
axion field could be defined. But that scenario would correspond to the decoupling 
of the axion from the SM sector which takes place in the massless quark case. 
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In conclusion, the axion coupling must be denned rigorously at the axion vacuum. 
However, above the TeV scale we can choose any basis for calculating the axino 
production if the contribution for E <1 TeV is not important. 

After the U(l)x rotation, the SU(2) singlets (quarks and leptons) carry Q' 
charges and hence the aYY coupling will not be absent: 

a Y C aY Y D g /r . „s 

~fI/N~ <2 ' 9) 

where B denotes here the field strength of the U(1) Y gauge boson and B denotes 
its dual. Since B^ = A^cosOw + Z^smOw, we obtain 

OL Y COS 2 9 W CaYY ^ ^ /n - n v 

— — aF em F em , (2.10) 

J a/ M 

where, F em is the field strength of the electromagnetic gauge boson and F em is its 
dual. This implies C a YY = C a -yy- 

In the DFSZ model with (d c , e)-unification, C a yy = 8/3. In the KSVZ model 
for e Q = 0, -1/3, and 2/3, C aYY = 0, 2/3 and 8/3, respectively [47]. Below the QCD 
chiral symmetry-breaking scale, C a77 and C a YY are reduced by 1.92. 

For our purpose the most important coupling will be that of axino-gaugino- 
gauge boson interactions which can be derived from Eq. (2.8) and written in a more 
conventional way as a dimension-5 term in the Lagrangian: 

Here and below B denotes the bino, the fermionic partner of the U(l)y gauge boson 
B and g stands for the gluino. 

One can also think of terms involving dimension-4 operators coming, e.g., from 
the effective superpotential where ^ is one of the MSSM matter (super)fields. 

However, axino production processes coming from such terms will be suppressed 
at high energies with respect to processes involving Eq. (2.11) by a factor m\/s, 
where s is the square of the center of mass energy. We will comment on the role of 
dimension-4 operators again below but, for the most part, mostly concentrate on the 
processes involving axino interactions with gauginos and gauge bosons, Eq. (2.11), 
which are both model- independent and dominant. 



3. Thermal production 

As stated in the Introduction, thermal production proceeds through collisions and 
decays of particles that still remain in the thermal-bath. As we will see, its efficiency 
will strongly depend on effective interaction strengths of the processes involved and 
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on characteristic temperatures of the bath. Particles like axinos and gravitinos are 
somewhat special, in the sense that their interactions with other particles are very 
strongly suppressed with respect to the Standard Model interaction strengths. There- 
fore such particles remain in thermal equilibrium only at very high temperatures. In 
the particular case of axinos (as well as axions and saxions), their initial thermal 
populations decouple at [13] 

At such high temperatures, the axino number density is the same as the one of 
photons and other relativistic species. In other words, such primordial axinos freeze 
out as relativistic particles. Rajagopal, Turner and Wilczek (RTW) [13] pointed out 
that, in the absence of a subsequent period of inflation, the requirement that the 
axino energy density be not too large (f^ < 1) leads to 

m 2 <12.8evf^)y (3.2) 
V. 9cs J 

where g c g = 1.5 and g* is the number of effectively massless degrees of freedom 
(particles with mass much smaller than temperature). In the MSSM at temperatures 
much higher than Msusy, one has g* = 915/4. The RTW bound (3.2) then takes a 
more well-known form [13] 

m s <2keV, (3.3) 

and the corresponding axinos would be light and would provide warm or even hot 
dark matter [13]. We will not consider this case in the following, primarily because 
we are interested in cold DM axinos. We will therefore assume that the initial 
population of axinos (and other relics, such as gravitinos), which were present in the 
early Universe, was subsequently diluted away by an intervening inflationary stage 
and that the reheating temperature after inflation was smaller than Td. It also had to 
be less than f a , otherwise the PQ symmetry would have been restored, thus leading 
to the well-known domain-wall problem associated with global symmetries. 

At lower temperatures, T < Tq, the Universe can be re-populated with axinos 
(and gravitinos) through scattering and decay processes involving superpartners in 
the thermal-bath. As long as the axino number density n„ is much smaller than 
n 7 , the number density of photons in thermal equilibrium, its time evolution will be 
adequately described by the Boltzmann equation: 

+ 3Hria = y](<r(z + j -> a H )u re i)n i n J - + J^(r(i -> a H ))rn. (3.4) 

Here H is the Hubble parameter, H(T) = a/ (tt 2 ^) / (90Mp) T 2 , where g* has been 
defined above and cr(i+j — > a+- ■ •) is the scattering cross section for particles i, j into 
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final states involving axinos, v re \ is their relative velocity, rii is the ith particle number 
density in the thermal-bath, T(i — > a+- ■ •) is the decay width of the ith particle, (• • •) 
stands for thermal averaging. (Averaging over initial spins and summing over final 
spins is understood.) Note that on the r.h.s. we have neglected inverse processes, 
since they are suppressed by n^. 

Solving the Boltzmann Equation. In order to solve the Boltzmann Eq. (3.4), 
it is convenient to introduce the axino TP yield 



TP 

Y~r = ^, (3.5) 



where s = (2ir 2 /45)g s *T 3 is the entropy density, and normally g s * = g* in the early 
Universe. We also change variables from the cosmic time t to the temperature T by 
—dt/dT = 1/HT for the radiation dominated era. One integrates the Boltzmann 
equation from the reheating temperature Tr after inflation down to zero. The yield 
can be written as 

^ TP = E F *:r t+ E r * dec ' ( 3 - 6 ) 

i,j i 

where the summation over indices i,j is the same as in Eq. (3.4). The expressions 
for Yj^ at and Yf cc are given by 



r 

^scat / 

i,3 ~ / 
JO 



R 



and 



>7-= r^ (r( '7/ r ' • •»>"■ ■ (3.8) 



o 



Explicit formulae for Yf^ 1 and Y^ cc can be found, for instance, in Ref. 

Scatterings. In the case of axinos the main productions channels are the scat- 
terings of (s)particles described by a dimension-5 axino-gaugino-gauge boson term 
in the Lagrangian (2.11). Because of the relative strength of a S) the most important 
contributions will come from strongly interacting processes. We will discuss them 
first. 

The relevant 2-body processes for strongly interacting particles i,j into several 
final states involving axinos are listed in Table 1. The corresponding cross sections 
<j n = a{i + j — > a + • • •), where the label n = A, . . . , J counts the different allowed 
combinations of the initial and final-state particles, can be written as 

a 3 

a n (s) = o<J n (s) (3.9) 

where y/s is the energy in the center-of-mass frame. Also listed there are the respec- 
tive: spin factor n spin (the number of spin combinations in the initial state), flavor 
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n 


Process 




^■spin 


np 


V1V2 


A 


g a + g b - a + g c 


!l /a6 c 


2 


4 


1 


1 


B 


g a + g b — > a + g c 


5 1 Xa6 
161 J 


'I 2 [log (slmlcr) - ^1 

1 O V / cn / 5 


4 


1 


3 


C 


g a + h -> a + <?j 


1 1 

8 1 


T a 1 


2 


2 


Np x 2 


1 


D 


5 ,a + <?fc a + <?j 


1 

32 


\ r pa 

\ I ,ik 


2 


4 


JVp x 2 


3 

/I 


E 


9j + <?fc -> a + g a 


1 

16 


\T a 
\ I jk 


2 


2 


iV F x 2 


3 


F 


g a + — > a + g c 


\\f abC 


2 [log (s/m 2 ff ) - f|] 


4 


1 


3 3 

4 4 


G 


g a + qk -> a + gj 


1 1 
4 1 


T a 1 
jk\ 


2 [log ( S /m 2 ff ) - 2] 


4 


iV F 


3 3 

4 4 


H 


<? a + <?fc — > a + 


1 1 
4 1 




2 [log ( S /m 2 ff ) - f ] 


2 


Np x 2 


3 
4 


I 


g fc + ^ -> a + g a 


1 

24 


\ r Pa 
l-Ljk 


2 


4 


Np 


3 3 

4 4 


J 


qk + qj a + <? a 


1 

24 


r-pa 
1 jk 


2 


1 


N F x2 


1 



Table 1: The cross sections for the different axino thermal-production channels involving 
strong interactions. Masses of particles have been neglected, except for the plasmon mass 
m e ff. See text for explanation of the different symbols. 

factor n F (number of color-triplet chiral multiplets), and number density factor rji, 
i — 1,2 (rji — 3/4 (1) for each initial-state fermion (boson). The group factors 
jabc an( j rpa^ Q £ g au g e g roU p SU(N) satisfy the usual relations ^ abc \f ahc \ 2 = 
N(N 2 - 1) and £ a £ . fe |T^| 2 = (iV 2 - l)/2. 

The diagrams listed in the Table are analogous to those involving gravitino pro- 
duction considered by Moroi et al. [8], and later by Bolz et al. [9], and we follow 
their classification. This analogy should not be surprising, since both particles are 
neutral Majorana superpartners. (We will return to the gravitino-axino analogy in 
Section 6.) In Fig. 1 we present the Feynman diagrams corresponding to channels A 
and G. The diagrams for the other channels can be drawn in a similar fashion. 

The diagrams in channels B, F, G and H involve t-channel exchange of massless 
gluons and are therefore divergent. In order to provide a mass regulator, one can 
introduce a plasmon mass m e fj representing the effective gluon mass due to plasma 
effect. We follow the prescription of Ref. [7], which has also been used in Refs. [8, 9], 
and assume m 2 ff = g 2 s T 2 where g s is the strong coupling constant. 

Note that we have neglected in the gluon thermal mass a potentially sizeable 
factor which counts the number of colored degrees of freedom of the plasma but kept 
a constant term in the logarithmically divergent cross sections, which depends on the 
procedure used to regulate the divergence. While a more proper treatment exists in 
the literature [49, 11], based on Thermal Field Theory and taking into account also 
the Fermi or Bose nature of the particles in the thermal-bath, we believe that the 
above approximation is good enough for our purposes. This belief is supported in 
particular by the fact that in the case of gravitinos, the full treatment agrees with 
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CHANNEL A 




Figure 1: Feynman diagrams contributing to channels A and G in Table 1. Thick dot 
denotes effective axino couplings to ordinary colored particles. 

an approximate one within a factor 3 [11]. For simplicity, in computing a n (s) we 
neglected all mass terms (involving gauginos, axinos and scalars) and kept only the 
plasmon mass ra^. 

Finally, we comment on the scattering processes involving the hypercharge mul- 
tiplet. It is clear that these will always be subdominant. This is not only a result of 
a weaker interaction strength relative to gluons and gluinos (so long as C a YY is not 
too large), but is also caused by a smaller number of production channels. Indeed, 
since U(l)y is Abelian, in this case many partial cross sections (all channel A, B and 
F contributions) vanish automatically. The others are given by the same expressions 
as in Table 1, with the substitution a\ — > oc y C^ Y y an d \T^\ 2 — > Yf, where Y{ is the 
hypercharge of the initial (and final) state. 
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Decays. In addition to scattering processes, axinos can also be produced through 
decays of heavier superpartners in thermal plasma. At temperatures T > rrig, these 
are dominated by the decays of gluinos into LSP axinos and gluons. The relevant 
decay width is given by 

and one should sum over the color indices a,b — 1, • • • , 8. 

In addition, at temperatures m x < T R < rrig, neutralino decays to axinos also 
contribute, while at higher temperatures they are sub-dominant. The relevant con- 
tribution is given by the decay of the bino component: 

2 r?2 ™3 / 9 




r(xi^a + B)= - a T l£ Z,2 U-ff- • (3-11) 



128tt3 ( /a /A0 



Here a em is the electromagnetic coupling strength, C aXi B = ZibC a YY / cos 2 6w where 
ZiB is the bino component of the ith neutralino x% (i = 1> 2, 3, 4). We use the basis 
Xi = ZnB + Z i2 W^ + Z i3 H® + Zi 4 H^ of the respective fermionic partners (denoted 
by a tilde) of the electrically neutral gauge bosons B and W3, and the MSSM Higgs 
bosons Hf, and H t . 

Discussion and Results. We have evaluated the integrals (3.7) and (3.8) numer- 
ically using expressions (3.9), (3.10) and (3.11). The results are presented in Fig. 2 
for representative values of f a = 10 11 GeV and rrig — rrig — 1 TeV. The respective 
contributions due to scattering as well as gluino and neutralino decays are marked 
by dashed, dash-dotted and dotted lines. 

It is clear that at high enough T R , much above rrig and rrig, scattering processes 
involving such particles dominate the axino production. For Tr ^> mg,m g , (a n v) is 
almost constant and thus F scat ~2x 10~ 5 MpT R x (%2 n &n)- 111 other words, F scat 
grows linearly as Tr becomes larger. (In the numerical calculation we have also taken 
into account the running of the strong coupling constant. Then y scat is given, to a 
very good approximation, by the same formula with a s replaced by q; s (Tr), and it 
grows like a 3 s (T R )T R . Using a s (T R ) instead of its value at M z gives a correction of 
up to a factor 10. For example, a s (M^) 3 /o; s (10 8 GeV) 3 = 5.75.) A similar result 
was estimated and used previously for the saxion in Ref. [26]. In contrast, the decay 
contribution above the gluino mass threshold, F dcc ~ 5 x 10~ 4 (M P r 5 /m|), remains 
independent of T R . This is not surprising since the scattering term in Eq. (3.4) 
behaves like const x riiiij ~ const x T 6 relative to the decay term which scales 
like m~rig ~ m~T 3 plus additional suppression (T) ~ T _1 . In the region of very high 
T R ~ 10 9 GeV, Y^ p becomes comparable with the axino yield in thermal equilibrium 
y^Q ~ 2 x 10~ 3 . Since, as mentioned below Eq. (3.4), we have neglected processes 
of axino re-annihilation to thermal-bath particles, there appears an artificial "edge" 
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Figure 2: Y~ TP as a function of Tr for representative values of f a = 10 11 GeV and m-q = 
m~g = 1 TeV. 

between the sloping solid curve representing Y~^ v and the dashed horizontal line 
corresponding to the axino equilibrium value ~ 2 x 10 -3 . Including axino re- 
annihilation would have the effect of rounding it. We leave it as is because we are 
not interested in very large values of T R . 

At T R roughly below the mass of the squarks and gluinos, their thermal popu- 
lation starts to become strongly suppressed by the Boltzmann factor e _m//T , hence 
causing a distinct knee in the scattering contribution in Fig. 2. It is in this re- 
gion that gluino decays (dash-dotted line) given by Eq. (3.10) become dominant, 
before they also become suppressed by the Boltzmann factor due to the gluino mass. 
For m x < Tr < mq,m g , the axino yield is well approximated by Y TP ~ Y dcc ~ 
5 x 10~ 4 (Mpr^/T^) e -m§ / TR , and depends sensitively on the reheating temperature. 

At still lower temperatures the population of strongly interacting sparticles be- 
comes so tiny that at Tr ~ m x neutralino decays given by Eq. (3.11) start playing 
some role, until they too become suppressed by the Boltzmann factor. We indicate 
this by plotting in Fig. 2 the contribution of the lightest neutralino (dotted line). By 
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comparing Eqs. (3.10) and (3.11) we can easily estimate that the bino-like neutralino 
contribution is suppressed by {a em /d>a^) 2 (m x /mg) 3 ' ~ 10~ 4 (with the color factor of 
8 included), assuming the usual gaugino mass relations of the MSSM. It is clear that 
the values of Y~f v in this region are so small that, as we will see later, they will play 
no role in further discussions. We therefore do not present the effect of the decay of 
the heavier neutralinos. 

As noted at the end of Section 2, there are also dimension-4 operators contribut- 
ing to axino production processes. One such operator is given by C a qq(rnq/ faja'jsqq* , 
where C aqq - is an effective axino-quark-squark coupling whose size is model-dependent, 
but which usually arises at the two-loop level. It contributes to processes involving 
quarks in Table 1. For example, in channel G, in addition to the t-channel gluon 
exchange diagram drawn in Fig. 1, there will now be a diagram with a squark ex- 
change in the s-channel. It will contribute a s {C aqq mq/ f a ) 2 {l/ s) compared to c? s j f 2 
of dimension-5 operators. Hence it is suppressed at high energies, s ^> m~, by a fac- 
tor {C^q/ a s) 2 {my s), and it can be non-negligible only around the mass threshold 
s ~ m|. This will, however, have a relatively small effect on the integration of the 
Boltzmann equation, unless Tr ~ m,q. For the DFSZ case, there is an additional 
contribution to the axino production coming from diagrams involving the Higgs su- 
permultiplet in the initial or final state; we neglect such part in the present study 
since it is model-dependent and we expect it to be subdominant. 

The sensitivity of Y~- p to the Peccei-Quinn scale f a is presented in Fig. 3. As is 
clear from Eq. (3.9), the rate of axino production is inversely proportional to f 2 . 

We emphasize that axinos produced in this way are already out of equilibrium. 
Their number density is very much smaller than n 7 (except T R ~ 10 9 GeV and above) 
and cross sections for axino re-annihilation into other particles are greatly suppressed. 
This is why in Eq. (3.4) we have neglected such processes. Nevertheless, even though 
axinos never reach equilibrium, their number density may be large enough to give 
Qa ~ 1 for large enough axino masses (keV to GeV range), as we will see later. 

Before closing this section, it is worth discussing how the presence of saxions 
could influence the above results. The saxions will be produced thermally in a way 
analogous to the axinos but next they will decay with a relatively short lifetime r s 
given by 

— - 6 -(l#iv^) 2 (l^)" 

This time needs to be compared with the time t s = 2.9 x 10~ 5 (1 TeV/m s ) 2 when the 
saxions dominate the energy density. If r s > t s , the saxion decay generates significant 
entropy and hence dilutes the particle species decoupled before the decay, such as 
the axions and the axinos. In other words, the present axion density is lowered 
and therefore the upper bound on f a is relaxed to values well beyond 10 12 GeV. 
Furthermore, the photon temperature now decreases more slowly, as pointed out 
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Figure 3: Y~- p as a function of Tr for representative values of f a and nig = nig = 1 TeV. 

in Ref. [26]. These effects become important once the saxion mass is greater than 
5 TeV [27] and f a is large. In this study, however, we will neglect the reheating 
effect due to saxion decay since we will assume that the saxion mass is below 1 TeV. 
Finally, in case the saxion mass is much smaller than the SUSY-breaking scale, as 
for instance in no-scale models, then the saxion decays very late and can play the 
role of a late-decaying particle [26] instead of generating entropy. 

4. Non-Thermal Production 

As discussed in the Introduction, axinos may also be produced in decay processes 
of particles which themselves are out of equilibrium, the decaying particle being one 
of the ordinary superpartners, the gravitino or the inflaton field. Below we will 
concentrate on the first possibility. 

Let one of the ordinary superpartners be the LOSP and the NLSP (next-to- 
lightest supersymmetric particle). (They do not have to be the same. The case 
when the role of the NLSP will be taken by the gravitino will be discussed below.) 
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A natural, albeit not unique, candidate for the LOSP is the lightest neutralino. 
This is because its mass is often well approximated by the bino mass parameter Mi 
(nearly pure bino case) or by the /x-parameter (higgsino limit), neither of which grows 
much when evolved from the unification scale down to the electroweak scale. It is 
therefore natural to expect the neutralino to be the LOSP. Furthermore, in models 
employing full unification of superpartner masses (such as the CMSSM/mSUGRA), 
a mechanism of electroweak symmetry breaking by radiative corrections typically 
implies fi 2 3> M\ . As a result, the bino-like neutralino often emerges as the lightest 
ordinary superpartner [46, 19]. 

Axino production from bino-like neutralino decay was analyzed in the previous 
paper [1]. For the sake of completeness, we summarize here the relevant results as well 
as make additional points regarding the case of a general type of neutralinos. The 
process involves two steps. The LOSPs first freeze out of thermal equilibrium and 
next decay into axinos. Since the LOSPs are no longer in equilibrium, the Boltzmann 
suppression factor in this case does not depend on the thermal-bath temperature and 
has to be evaluated at the LOSP freeze-out temperature, which is well approximated 
by T f ~ m x /20. At T < T f and for Y X <^H but large enough for the decay to occur 
in the radiation-dominated era, one finds that Y x is roughly given by [1] 

" Tf dVm 2 (T x ) T , 



r x (T)^y E Q(T f )ex P 



T'3 H(m x ) 



(4.1) 



where y^ Q (T f ) contains the Boltzmann suppression factor evaluated at T f , (T X ) T is 
the thermally averaged decay rate for the neutralino at temperature T, and H(m x ) 
is given below Eq. (3.4). While we have used here the neutralino as the LOSP, the 
same mechanism would work for other ordinary superpartners as well. 

It is clear that, in order for a two-step process to occur, the decay width of 
the LOSP has to be sufficiently small to allow for the freeze-out in the first place. 
Fortunately this is what usually happens since the interactions between the LOSP 
and the axino are suppressed by the large scale f a . On the other hand, the lifetime 
of the LOSP must not be too large, otherwise the decay into axinos and ordinary 
particles would take place too late, during or after nucleosynthesis, and could destroy 
successful predictions for the abundance of light elements. This will be discussed in 
more detail in the next section. Here we note that it is truly remarkable that the 
relative strength of the axino interaction, in comparison with SM interactions, is 
such that the decay width falls naturally between these two limits. 

It is worth mentioning that in principle the LOSP may decay into axinos even 
in the absence of freeze-out. It has also been recently pointed out in Ref. [50] that 
even in the case of very low reheating temperatures Tr, below the LOSP freeze-out 
temperature, a significant population of them will be generated during the reheating 
phase. Such LOSPs would then also decay into axinos as above. The process will 
be non-thermal in the sense that the decaying LOSPs will not be (yet) in thermal 
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equilibrium. We will not pursue this possibility here but comment on it again when 
we discuss our results. 

In Ref. [1] we considered the non-equilibrium process 

X^a 7 . (4.2) 

This decay channel is always allowed. The decay rate for process (4.2) can be easily 
derived from Eq. (3.11): 



128jt" (fJNY \ m 



where C axi = (C a yy / cos 9w)Zn, with Zu standing for the bino part of the lightest 
neutralino 4 . 

The neutralino lifetime can be written as 

^^-^c^zfXTm) lira) {—) i 1 -^) 

(4.4) 

For large enough neutralino masses, an additional decay channel into axino and 
Z opens up, 

r (x - aZ) = tan 2 9 w -^L- 2 x PS {% ^) , (4.5) 

1287r3 (fa/N) 2 \ml ml) 



where the phase-space factor is given by 



PS(x,y) = v / l-2(x + y) + (x-y) 



x 2 



:i-yY--(i-QVy + y)-Y 



(4.6) 



Note that this channel is always subdominant relative to x ~^ ^7 because of 
both the phase-space suppression and the additional factor of tan 2 9w- As a result, 
even at m x ^> mz,rria, t(x — * aZ) ~ 3.35 t(x — > aj). It is also clear that the 
neutralino lifetime rapidly decreases with its mass (~ l/m x ). On the other hand, if 
the neutralino is not mostly a bino, its decay will be suppressed by the Z n -factor in 
C axi . 

Other decay channels are the decay into axino and Standard Model fermion pairs 
through virtual photon or Z, but they are negligible with respect to the previous 
ones. We will discuss them later since, for a low neutralino mass, i.e. long lifetime, 
they can, even if subdominant, produce dangerous hadronic showers during and after 
nucleosynthesis. 

In the DFSZ type of models, there exists an additional Higgs-higgsino-axino cou- 
pling, usually related to the MSSM //-term in simplest realizations. For example, by 



4 In Eq. (4.3) we have corrected an overall numerical factor with respect to Eq. (6) of Ref. [1]. 
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using the fields Si and £2, which break the PQ symmetry in the superpotential given 
by (2.5), one can generate the MSSM Higgs mass term through the renormalizable 
interactions 

W„ = X t S t H u H d (4.7) 

for i — 1, 2, with a very small A« ~ raw/ fa or through a non-renormalizable interac- 
tion like 

^ = Ai«. (4.8) 

After the PQ symmetry breaking, such superpotentials not only generate the 
MSSM //-term [35, 51], but in general also give rise to a Higgs-higgsino-axino cou- 
pling of order fi/f a - In the non-renormalizable case a four-fermion coupling of order 
/V/a i s generated. This coupling can be neglected since it does not contribute to 
the neutralino decay but the other can if the higgsino component of the neutralino is 
non- negligible. Note that after electroweak symmetry breaking such couplings pro- 
duce a mixing between the axino and the higgsinos but in general these are of the 
order of fj,v/ f a , and therefore much less than //, so that we can continue to consider 
the axino as an approximate mass eigenstate which is nearly decoupled from the 
MSSM. 

An example of such an effective MSSM+axino model has been recently analyzed 
in Ref. [52] for an effective potential of the type 

Wf = n (l + 1$) H u H d , (4.9) 

where e ~ v/ f a ~ 1CT 8 and v is the Higgs VEV of the order of the weak scale. One 
needs to perform a full diagonalization of the 5x5 neutralino+axino mass matrix. 
Since in this case the higgsinos couple more strongly (or rather, less weakly) to the 
axino, the non-negligible mixing between axino and higgsino introduces additional 
decay channels, e.g. the one with intermediate virtual right-handed sleptons, that 
are absent in the KSVZ case. 

By rescaling the results of Fig. 2 of [52] to our central value f a = 10 11 GeV, one 
finds that for low bino mass, m x < 120 GeV, = 50 GeV, and specific choices of 
other supersymmetric parameters, in particular masses of the right-handed sleptons 
degenerate and larger than m x , the decay time of the neutralino into axino and 
lepton or quarks pairs is of the order of 

t(x — > all and aqq) ~ 0.02 sec; (4-10) 

it decreases at higher masses because of the dependence on m x and of the opening 
of new channels (bb and WW*), down to 3 x 10~ 6 sec at m x = 200 GeV. We then see 
that even in the DFSZ case the lifetime is of the right order of magnitude to allow 
the freeze-out and decay process since H(m x ) ~ 2 x 10 7 sec _1 (m x /100 GeV) 2 . 
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Of course this result strongly depends on the type of DFSZ model considered: 
for example, the coupling e between the axino and the Higgses can be suppressed or 
enhanced by mixing angles in the PQ sector since the axion multiplet is in general 
a combination of different multiplets (Si and S 2 in the above example, compare 
Eq. (2.5)). Such angles are assumed to be of order 1 in the simple estimate above. 

We therefore conclude that DFSZ-type models have to be analyzed case by case, 
but in general we expect them to give an even more efficient implementation of the 
non-thermal production through LOSP decay. A scenario of this type has been re- 
cently studied in detail in [40] with the conclusion that an axino (or better in this case 
flatino) LSP and CDM candidate is still possible and that its population is produced 
naturally after a phase of thermal inflation, with very low reheat temperature. 

While in the discussion of the non-thermal production of axinos, we have con- 
centrated on the neutralino as, in some sense, the most natural choice for a parent 
LOSP, we reiterate that in principle one could also consider other choices for the 
LOSP, including charged particles, which we will not do here. 

Finally, we note for completeness that another way of producing axinos non- 
thermally is to consider gravitino decays [25]. In this case it is the gravitino, and not 
the LOSP, that is the NLSP. Such axinos are produced very late, at 10 8 sec. As we 
have stated in the Introduction, this channel may solve the gravitino problem [25]. 
We will come back to this point in the next Section. Furthermore, axinos may be 
produced during (p)reheating in the decay of the inflaton field. Such processes are 
strongly model-dependent and we will not consider them here. 



5. Constraints 

Several non-trivial conditions have to be satisfied in order for axinos to be a viable 
CDM candidate. First, we will expect their relic abundance to be large enough, 
Qah 2 — 0.2. This obvious condition will have a strong impact on other bounds. 
Next, the axinos generated through both TP and NTP will in most cases be initially 
relativistic. We will therefore require that they become non-relativistic, or cold, much 
before the era of matter dominance. Furthermore, since NTP axinos will be produced 
near the time of BBN, we will require that they do not contribute too much relativistic 
energy density to radiation during BBN. Finally, associated decay products of axino 
production will often result in electromagnetic and hadronic showers which, if too 
large, would cause too much destruction of light elements. In deriving all of these 
conditions, except for the first one, the lifetime of the parent LOSP will be of crucial 
importance. 

First, from now on we will assume that the axinos give a dominant contribution 
to the matter density at the present time. This can be expressed as 

Skh 2 = m~ a Y~ a ~ 0.2, (5.1) 
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where s is given under Eq. (3.5) and we have assumed no significant entropy produc- 
tion. Since p crit /h 2 = 0.8 x 10~ 46 GeV 4 and </ s *(T now ) = 3.91, this can be re-expressed 
as 

m~ a Y~ a ~ 0.72 eV ( -^-\ , (5.2) 

which readily applies to both TP and NTP relics. 

We note in passing that, for the initial population of axinos, the yield at decou- 
pling is approximately Ya — 5^ ~ 2 x 10~ 3 , which gives 

m s ~0.36keVf-^-J . (5.3) 

This is an updated value for the Rajagopal-Turner-Wilczek bound (3.3). 

For a fixed value of fl^h 2 , lower values of Ya give correspondingly larger m^. For 
example, the turn-over region of Y& ~ 10~ 10 in Fig. 2, where gluino decay becomes 
more important than scattering, gives m„ ~ 7.2 GeV (Qah 2 /0.2). 

Cold Axinos Next, we want to determine the temperature of the Universe at 
which the axinos will become non-relativistic. In nearly all cases axinos are initially 
relativistic and, due to expansion, will become non-relativistic at some later epoch, 
which depends on their mass and production mechanism. As mentioned in the Intro- 
duction, in the case of thermal production, even though the axinos are not in thermal 
equilibrium, they are produced in kinetic equilibrium with the thermal-bath. Their 
momenta therefore have a thermal spectrum, derived from the scattering particles in 
the plasma. The axinos become non-relativistic when the thermal-bath temperature 
reaches the axino mass: 

T N r ^ ma, (5.4) 

even though they are not in thermal equilibrium. 

In the case of non-thermal production, the momentum of the axinos depends 
strongly on the production mechanism. Here we will consider, as before, only the 
production through out-of-equilibrium neutralinos. In that case, axinos will be pro- 
duced basically monochromatically, all with the same energy, roughly given by m x /2, 
unless they are nearly mass-degenerate with the neutralinos. This is so because the 
neutralinos, when they decay, are themselves already non-relativistic. Thus the ax- 
inos will normally become non-relativistic only at later times, through momentum 
red-shift. However, in this case the simple estimate (5.4) used above for Tnr does 
not hold, since the axinos are not in kinetic equilibrium with the thermal-bath. 

First we need to determine the temperature T^ cc corresponding to the time of 
neutralino decay tdcc- We will make the approximation that most axinos are produced 
in neutralino decays at the time equal to the neutralino lifetime r x = l/r x , tdcc — t x 
(a sudden-decay approximation), which is a very reasonable one. 

Because not all neutralinos decay at the same time as assumed in the sudden- 
decay approximation, at a given time the momenta of produced axinos will actually 
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have a distribution that will be limited from above by m x /2. Let us find this momen- 
tum distribution fa{to,p) of axinos coming from neutralino decays. When the axino 
was produced by the decay of neutralino at time t, it had the momentum p = m x /2, 
and would be next red-shifted while the axino was relativistic. Since the scale factor 
R grows as R oc t 1 ^ 2 during the RD epoch, at a later time to, the axino produced at 
a time t would have the momentum 



KMo) 



TO. 



( R(t) \ 



TO, 



1/2 



(5.5) 



2 \R(t Q )J 2 

By counting the number of produced axinos during the time interval dt at t, we get 



N x — = fa(t ,p)dp, 



(5.6) 



where N x (t) = N x0 e t ^ Tx is the number of neutralinos at t which is obtained from 
the decay equation dN x /dt = (1/t x )N x . Thus, 

2p 



. N x dt 
t x dp 



iV v 



From this, we can find the total energy of axinos 



-( P /p(T x ,to)Y 



(5.7) 
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(5.8) 



For to ^> t x , the energy density of axinos becomes 



Pa = 



Ea(t ) _ Vrr m x f t. 
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This differs from the sudden-decay approximation only by a factor ps 0.89 

Using the sudden-decay approximation we obtain 
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(5.10) 
(5.11) 



where we have taken (?*(T ( j ec ) ~ 10. 

Axino momenta will red-shift with temperature as 



p(T) 



TO^ 
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(5.12) 



as long as the axinos remain relativistic. (The factor involving g Si , at different times 
accounts for a possible reheating of the thermal-bath when some species become 
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non-relativistic.) Axinos become non-relativistic at the epoch when p(Tnr) — ma, 
which gives 



NR 



m x \9s*( T nr) 



1 V; :r , (5.13) 



100GeV\ / T dcc 



= 2 - 7xi °-M^f (514) 

where we have taken g s *(T dcc ) ~ g*(T dec ) ~ 10 and g s *(T NR ) = 3.91. 
Using an explicit expression for T dcc (Eq. (5.11)) then leads to 

^♦(^dec)^ 1//3 CaYF Zu ( m a 



= 53.8 keV 



g s *(T NR )J gl/\ Tdcc ) 



llGeW UoGeW V /«/# 7 j 



For T dcc < 200 MeV, <? s *(T dec ) ~ ^(T dcc ) ~ 10 and for T < 500 keV, (? S *(T) ~ 3.91, 
which finally leads to 

T NR = 4.2 x (^) " 2 (^) ■ (5.17, 

This epoch has to be compared with the matter-radiation equality epoch. As- 
suming that axinos constitute the largest part of the matter density, Pmatter(T'cq) = 
ps(T cq ), and since 

27T 2 

p~ a = m~ a Yas(T) = m ~Ya—g s *(T) T 3 , (5.18) 
the temperature at matter-radiation equality is given by 

Te, = (5.19) 

^g*\-L cqj 

Using Eq. (5.2) and assuming g s *(T cq ) = 3.91 and g*(T eq ) = 3.36 allows us to express 
Eq. (5.19) as 



T eq = l.leV(^-^-j, (5.20) 

which holds for both thermal and non-thermal production. 

In the former case, by comparing with Eq. (5.4), one can easily see that T NR > T cq 
is satisfied for any interesting range of mj. In the case of NTP, in Eq. (5.19) we 
make a substitution i~ = Y X (T { ), the neutralino yield at freeze-out. The condition 
Tnr 3> T eq is satisfied for 
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where the Td ec dependence on the neutralino lifetime is given in Eq. (5.11). If axinos 
were lighter than the bound (5.21), then the point of radiation-matter equality would 
be shifted to a later time around T NR . Note that in this case the axino would not 
constitute cold, but warm or hot dark matter. We will see, however, that in the NTP 
case discussed here other constraints would require the axino mass to be larger than 
the above bound, so that we can discard this possibility. 

BBN constraint on NTP axinos In the case of non-thermal production, most 
axinos will be produced only shortly before nucleosynthesis and, being still relativis- 
tic, may dump too much to the energy density during the formation of light elements. 
The axino energy density now reads 



f>a(T) = YaS(T) P (T) 

'2tt 2 \ / ^(T) 4 / 3 \ (m x Y z 



45 J VM7dec) 1/3 y V2T dcc 



V^(^dcc) i/3 y \m~a) v t Acc J v 0.2 ; 1 ; 

where Eqs. (5.2), (5.12), and an expression for the entropy density s given under 
Eq. (3.5) have been used. In order not to affect the Universe's expansion during 
BBN, the axino contribution to the energy density should satisfy 

^ < 5N U , (5.24) 

Pv 

where p u = (7r 2 /30)(7/4)T 4 is the energy density of one neutrino species. Agreement 
with observations of light elements requires [53] SN U = 0.2 to 1. Using Eq. (5.11) 
leads, after some simple algebra, to 

^ nrrA , ^ 1 / m v \ (lMeV\ fttah 2 
m~ n > 274 keV 



llOOGevJ [-T— (5 - 25) 



SN V VlOOGeV/ V. T dec J \ 0.2 
. 181lfV 1 1 / 100GeV \ 1/2 / f a /N \ fn~ a h 2 \ 

where we have chosen the most conservative case g s *(T) = g s *(T dec ) ~ 10. (For 
g s *(T) = 3.91 the bounds become respectively 78keV and 55keV.) As before, Td ec 
depends on the neutralino lifetime via Eq. (5.11). 

Showers In the NTP case, if neutralino decays take place during or after BBN, 
produced bosons may lead to a significant depletion of primordial elements [18]. One 
often applies the crude constraint that the lifetime be less than about 1 sec. This 
in our case would provide a lower bound on m x . A detailed analysis [53] provides 
limits on the abundance of the decaying particle versus its lifetime. (See, e.g., Fig. 3 
of Ref. [53].) 



28 



First, the photons produced in reaction (4.2) carry a large amount of energy, 
roughly m x /2. If the decay takes place before BBN, the photon will rapidly ther- 
malize via multiple scatterings from background electrons and positrons (7 + e — > 
7 + 7 + e) [54, 15]. The process will be particularly efficient at plasma temperatures 
above 1 MeV, which is the threshold for background ee pair annihilation, and which, 
incidentally, coincides with a time of about 1 sec. But a closer examination shows 
that also scattering with the high-energy tail of the CMBR thermalize photons very 
efficiently. As a result, the decay lifetime into photons can be as large as 10 4 sec. By 
comparing this with Eq. (4.4) we find that, in the gaugino regime, this can be easily 
satisfied for m x < mz- It is only in a nearly pure higgsino case and for a mass of 
tens of GeV that the bound would become constraining. We are not interested in 
such light higgsinos for other reasons, as will be explained later. 

A much more stringent constraint comes from considering hadronic showers from 
gg-pairs. These will be produced through a virtual photon and Z exchange, and, 
above the kinematic threshold for x ~^ also through the exchange of a real 
Z-boson. We will now discuss this in some detail. 

By comparing with Fig. 3 of Ref. [53], we can see that the bound on hadronic 
showers can be written as 



xBR( X ^qq)<F(r x ), (5.27) 



where BR(x —> qq) is the branching ratio of neutralinos into axinos plus gg-pairs 
and F{t x ) can be read out from Fig. 3 of Ref. [53]. 
Since 1^ — ^x(Tf), an< ^ remembering that 



n x _ 71 9s* 



F X ~7.04Y X , (5.28) 



n 7 45C(3) 
one can express the condition (5.27) as 

As stated above, one can write 

BR(x — > gg) = BR(x — > 7* — > qq) + BR(x — ► Z/Z* — > gg) + interference term. 

(5.30) 

Assuming <C m x ,mz, in the region m x < mz, we find that BR(x —> 7* —> 
qq) dominates and (summing over all light quark species up to the bottom) is of the 
order of 0.03-0.04, slowly increasing with m x . Above the Z threshold the intermedi- 
ate Z channel also becomes sizeable, and the branching ratio grows faster with m x , 
up to ~ 0.06 at m x = 150 GeV. The interference part is always negligible. 

We have computed BR(x — * qq) both below and above the threshold for the 
process x ~^ f° r negligible axino mass and used it in our numerical analysis. Our 
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Figure 4: Lower bound on the axino mass from considering hadronic showers according 
to the condition (5.29), for C a YyZ\\ = 1 and f a /N = 10 11 GeV. The bound disappears for 
m x = 150 GeV when the lifetime drops below 0.1 sec. 

results are presented in Fig. 4 in the bino case Z\\C a YY — 1 and for f a /N = 10 11 GeV. 
For the general case, one can still read out the lower bound on the ratio rria/m x from 
Fig. 4 by replacing m x with (G (l yyZ 11 ) 2 / 3 (10 11 GeV / f a /N) 2 ^ 3 m x . In fact the main 
dependence on m x in Eq. (5.29) comes from the factor m x and the neutralino lifetime 
in F(t x ), while BR(x —> QQ) only slowly increases with m x . 

These bounds are clearly much more stringent than those in (5.21) and (5.25), but 
they are, at the same time, strongly sensitive to the neutralino mass and composition. 
Notice that, in the region of low m x , the bound on m„ gradually strengthens with 
increasing m x (and therefore decreasing neutralino lifetime) because, for lifetimes of 
order 2-20 sec, the function F(t x ) is shallow and the branching ratio and x both 
increase. However, as the lifetime drops below 1 sec, F(t x ) increases steeply and 
so the rria bound decreases almost linearly in m x before disappearing altogether for 
m x > 150 GeV (in the bino case) when r x < 0.1 sec. 

If the higgsino component of the decaying neutralino increases, so does the life- 
time. There are two points to note here. One is that the neutralino yield now 
becomes much smaller, owing to co-annihilation with the next to lightest neutralino 
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and lightest chargino until m x > 500 GeV, when the lifetime will be suppressed 
again. Furthermore, in the DFSZ model, new channels are present, as we discussed 
in the previous section, for which a typical lifetime will be very much smaller than 
lsec [52]. Only in the KSVZ model in the higgsino-like case, in a relatively small 
region 60 < m x < 150 GeV (the lower number being the current rough experimental 
bound), will there exist some restriction on the combination of neutralino mass and 
higgsino purity, which would further depend on the neutralino yield at freeze-out. 
Thus we conclude that overall we find no significant restriction on the neutralino 
mass from the lifetime constraint. 

In summary, a lower bound m„ > (9(300 keV) arises from requiring either the 
axinos to be cold at the time of matter dominance or that they do not contribute too 
much to the relativistic energy density during BBN. The constraint from hadronic 
destruction of light elements can be as strong as > 360 MeV (in the relatively 
light bino case), but it is highly model-dependent and disappears for larger m x . 

6. Thermal vs. Non-Thermal Production 

We now proceed to comparing the relic abundance of axinos produced non-thermally 
in neutralino decays with the thermal production case analyzed in Section 3. Clearly, 
in the TP case the axino yield is primarily determined by the reheating temperature 
(for a fixed f a ). For large enough Tr (Tr ^> m-g,mq), it is proportional to Tr//^. 
(Compare Figs. 2 and 3.) For the reheating temperature below the mass threshold 
of strongly interacting sparticles, it becomes Boltzmann suppressed by the factor 
~ exp- m 5/ TR (or ~ exp- m ?/ TR ). 

If the axino is the LSP, the present fraction of the axino energy density to the 
critical density is given by Eq. (5.2). This relation allows us to redisplay the results 
of Fig. 2 in the plane (m^T^). For thermally regenerated axinos, Eq. (5.2) gives an 
upper bound on Tr as a function of the axino mass from the requirement &~ P h 2 < 1. 
This bound is depicted by a thick solid line in Fig. 5. 

A digression is in order here. An expert reader will have noticed that, while in the 
axino case, the bound flah 2 < 1 gives Tr ~ 1/rria, in the gravitino case the analogous 
bound gives T R ~ mg. (See, e.g., Fig. 1 in Ref. [8].) This difference is caused by the 
fact that the crucial effective Lagrangian dimension-5 operator, which is responsible 
for the bound, is proportional to in the gravitino case, but exhibits no 

dependence in the case studied here. (Compare Eq. (2) or Ref. [8] with Eq. (2.8).) 

Since axinos couple like l/f a , one would naively expect that Tr for which TP 
becomes important (or when too much relic abundance is generated) would be just 
fa/Mp of that for gravitinos and thus hopelessly low, Tr < 10" 7 GeV. This is, 
however, not the case, since the gravitino production is dominated by the goldstino 
component, whose interaction is suppressed by the supersymmetry-breaking scale 
M s , rather than the Planck scale: for example, the coupling to the gluino is oc 
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Figure 5: The solid line gives the upper bound from thermal production on the reheating 
temperature as a function of the axino mass. The dark region is the region where non- 
thermal production can give cosmologically interesting results (f2~ TP /i 2 ~ 1) as explained 
in the text. We assume a bino-like neutralino with m x = 100 GeV and f a = 10 11 GeV. The 
region of Tr > Tf is somewhat uncertain and is shown with light-grey color. A sizeable 
abundance of neutralinos (and therefore axinos) is expected also for Tr < Tf [50] but 
has not been calculated. The vertical light-grey band indicates that a low range of mj 
corresponds to allowing SM superpartner masses in the multi-TeV range, as discussed in 
the text. The division of hot, warm and cold dark matter as a function of the axino mass 
shown in the lower left part is for axinos from non-thermal production. 

rrig/Mg. Thus T R at which TP becomes significant will be of order 10 3 -10 4 GeV, 
which is on the low side but still acceptable. 

In the NTP case, the yield of axinos is just the same as that of the decaying 
neutralinos. This leads to the following simple relation [1] 

Q~h 2 — — n x h 2 . (6.1) 
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We remind the reader that Q x h 2 stands for the abundance that the neutralinos 
would have had today, had they not decayed into axinos. It is related to the neu- 
tralino yield through an analogue of Eq. (5.2) (which actually applies to any stable 
massive species) and is determined by the effective cross section of neutralino pair- 
annihilation (as well as co-annihilation) into ordinary particles. When the rate of this 
process becomes less than the rate of the expansion of the Universe, the neutralinos 
freeze out. Typically this happens at freeze-out temperatures of T f ~ m x /20. 

In contrast to the TP case, the NTP axino yield will also be independent of 
the reheating temperature (so long as Tr ^> Tf). In order to be able to compare 
the two production mechanisms, we will therefore fix the neutralino mass at some 
typical value. Furthermore we will map out a cosmologically interesting range of 
axino masses for which fi~ TP ~ 1. 

Our results are presented in Fig. 5 in the case of a nearly pure bino. We also 
fix m x = 100 GeV and f a = 10 11 GeV. The dark region is derived in the following 
way. It is well known that Q x h 2 , the relic abundance of neutralinos, can take a 
wide range of values spanning several orders of magnitude. In the framework of the 
MSSM, which we have adopted, global scans give Q x h 2 < 10 4 in the bino region at 
m x ~ 100 GeV. (This limit decreases roughly linearly (on a log-log scale) down to 
~ 10 3 at m x ~ 400 GeV.) For m x = 100 GeV, by using Eq. (6.1) we find that the 
expectation f2~ TP /i 2 ~ 1 gives 

10MeV<m a <m x . (6.2) 

We note, however, that the upper bound Q x h 2 < 10 4 comes from allowing very 
large M S usy (i-e. sfermion and heavy Higgs masses) in the range of tens of TeV. 
Restricting all SUSY mass parameter below about 1 TeV reduces Q, x h 2 below 10 2 
and, accordingly, increases the lower bound mj > 1 GeV. For the sake of generality, 
in Fig. 5 we have kept the much more generous bound (6.2) but we marked a low 
range of with a light grey band to indicate the above point. 

Likewise, for reheating temperatures just above T f , standard estimates of Q x h 2 
become questionable. We have therefore indicated this range of Tr with again light 
grey color. It has also been recently pointed out in Ref. [50] that a significant 
population of LOSPs will be generated during the reheating phase even at T R below 
the LOSP freeze-out temperature. Such LOSPs would then also decay into axinos 
as above. We have not considered such cases in our analysis and accordingly left the 
region Tr < Tf blank, even though in principle we would expect some sizeable range 
of flah 2 there. 

We can see that, for large Tr, the TP mechanism is more important than the 
NTP one, as expected. Note also that in the TP case the cosmologically favored 
region (0.2 < Qah 2 < 0.4) would form a very narrow strip (not indicated in Fig. 5) 
just below the f2~ p = 1 boundary. In contrast, the NTP mechanism can give the 
cosmologically interesting range of axino's relic abundance for a relatively wide range 
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of rria so long as T R < 5 x 10 4 GeV. Perhaps in this sense, the NTP mechanism can 
be considered as somewhat more robust. 

We have also marked in Fig. 5 some of the bounds discussed in Section 5. They 
are normally not as restrictive as the shaded region in the Figure. For example, the 
ranges of where non-thermally produced axinos would be hot /warm/cold dark 
matter are denoted, following Eqs. (3.3) and (5.21). The potentially most stringent 
bound from hadronic showers would require ms > 360 MeV, but it is very sensitive 
to m x and disappears for m x > 150 GeV (compare Fig. 4). 

At larger bino mass, the lower bound on also increases, following Eq. (6.1) 
and because Q, x h 2 decreases with m x . Other bounds do not give any additional 
constraints. In the higgsino case, one typically has Q, x h 2 <C 1 (or, more properly, the 
higgsino number density at freeze-out is very much smaller than that of the bino with 
the same mass) owing mostly to co-annihilation, as mentioned at the end of the last 
section. A heavy higgsino with m x > 500 GeV nevertheless remains an option for 
a LOSP. In this case, however, additional model-dependent (dimension-4) operators 
will contribute to the TP mechanism, and the upper curve in Fig. 5 will probably be 
moved upwards. Overall we find the bino case to be much more natural and robust. 

7. Implications and Conclusions 

The intriguing possibility that the axino is the LSP and the dark matter WIMP 
possesses a number of very distinct features. This makes this case very different 
from those of both the neutralino and the gravitino. In particular, the axino can 
be a cold DM WIMP for a rather wide range of masses in the MeV to GeV range 
and for relatively low reheating temperatures T R < 5 x 10 4 GeV. As T R increases, 
thermal production of axinos starts dominating over non-thermal production and the 
axino typically becomes a warm DM relic with a mass broadly in the keV range. In 
contrast, the neutralino is typically a cold DM WIMP (although see Ref. [55]). 

Low reheating temperatures would favor baryogenesis at the electroweak scale. It 
would also alleviate the nagging "gravitino problem" . If additionally it is the axino 
that is the LSP and the gravitino is the NLSP, the gravitino problem is resolved 
altogether for both low and high T R . 

Phenomenologically, one faces a well-justified possibility that the bound Q x h 2 < 
1, which is often imposed in constraining a SUSY parameter space, may be readily 
avoided. In fact, the range Q x h 2 ^> 1 (and with it typically large masses of super- 
partners) would now be favored if the axino is to be a dominant component of DM 
in the Universe. Furthermore, the lightest ordinary superpartner could either be 
neutral or charged but would appear stable in collider searches. 

The axino, with its exceedingly tiny coupling to other matter, will be a real 
challenge to experimentalists. It is much more plausible that a supersymmetric 
particle and the axion will be found first. Unless the neutralino (or some other 
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WIMP) is detected in DM searches, the axino will remain an attractive and robust 
candidate for solving the outstanding puzzle of the nature of dark matter in the 
Universe. 
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